BELLE Data on the 7r 7*7 Form Factor: A Game Changer? 
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We extend our analysis [l(] of the 7r°7*7 form factor by including a comparison with the new 
BELLE data @]. The necessity of new precision measurements in a broad interval of momentum 
transfers is emphasized. 
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In this note we analyze the impact of the new BELLE 
data on the theoretical status of the 7r°7*7 form fac- 
tor in QCD and the constraints on the pion distribution 
amplitude (DA). It has been a hot subject over the last 
several years, with a strong scaling violation suggested 
by the BaBar measurements [3] fuelling a flurry of theo- 
retical activity. 

This analysis has to be viewed as an addendum to 
our work Ref. |l| where we have given an updated NLO 
analysis of the form factor in the framework of QCD 
collinear factorization, taking into account higher-twist 
corrections up to twist-six and also soft nonperturbativc 
corrections estimated using dispersion relations and du- 
ality. In this note we follow Ref. [l| closely both for the 
formalism and notations; We also refer to [1] for the com- 
plete list of relevant references. 

A compilation of the experimental data by BELLE, 
BaBar and CLEO collaborations 0-0] is shown in Fig. Q] 
together with QCD calculations using several models 
of the pion DA: asymptotic <j>^ s {x) — 6x(l — x) (solid 
line) , BMS model |5j (short dashes) , "holographic" model 
4% ol (x) = (8/tt) y/6x(l - x) U (long dashes), a model II 
of Ref. Q (dash-dotted) and "flat" ^(x) = 1 |7j (dots). 
One sees that the BELLE data are systematically lower 
than those of BaBar in a broad range of Q 2 and can be 
described e.g. by the "holographic" model. The quality 
of the fit is x 2 — 0-64 (per data point) adding theoretical 
uncertainties and experimental errors in quadrature. 

One can try to quantify the difference in the pion DAs 
which are required to describe these two data sets. To 
this end we take the model II of Ref. [l[ as an example 
(dash-dotted curve in Fig.[T]), and modify the parameters 
in order to obtain a good fit for BELLE data instead 
of BaBar. We end up with the following values of the 
Gegenbauer coefficients at the scale 1 GeV: 



a 2 = 0.10(0.14) , a 4 
a 6 = 0.10(0.18) , a 8 



0.10(0.23) , 
0.034(0.05), (1) 



where the numbers in parenthesis correspond to the 
model II of Ref. pj and, by design, describe the BaBar 
data Q . The transition form factor calculated using the 
pion DA defined by Eq. ([1]) is shown in Fig. [2] together 
with the estimated theoretical uncertainty. It is very sim- 
ilar to the "holographic" model, apart from being some- 
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FIG. 1: The pion transition form factor for the "asymp- 
totic" (solid line), "BMS"[fJ (short dashes), "holographic" 
(long dashes), "model II" of Ref . Q (dash-dotted) and "flat" 
(dots) pion DA. The experimental data are from (squares) , 
3] (circles) and [J] (open triangles). 
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FIG. 2: The pion transition form factor for the model of 
the pion DA described in the text. The estimated theoretical 
uncertainty is shown by the shaded area. The experimental 
data are from (squares) and [|| (open triangles). 



what lower at small Q 2 and provides a better overall fit in 
this region. Adding the theoretical uncertainty and the 
experimental errors in quadrature, we obtain for our new 
model x 2 — 0-44 (per BELLE data point). For compar- 
ison, we get x 2 — 1-06 for the description of the BaBar 
data in the model II of Ref. [l|. These two models for 
the pion DA, describing alternatively either BELLE or 
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2. intermediate Q 2 ~ 10 GeV 2 




FIG. 3: The model of the pion DA at the scale 1 GeV used 
in Fig. [2] (solid line) compared to the model II of Ref. [l[ 
shown by the dash-dotted curve in Fig. [TJ (dashed) and the 
asymptotic DA (dotted). 



BaBar data, are compared with each other in Fig. [3] It 
is seen that the BaBar measurements imply an enhance- 
ment of the DA in the end-point region, but apart from 
that the differences are minor. This is not surprising since 
the overall discrepancy between the BELLE and BaBar 
data sets is within 1.5-2 standard deviations. 

We conclude that more precise form factor measure- 
ments for both small and large momentum transfers are 
needed in order to arrive at a definite conclusion. Let us 
discuss the role of different regions. 



1. moderate 



2-6 GeV 2 



The QCD predictions in the Q 2 ~ 2 - 6 GeV 2 rely 
rather heavily on the duality assumption for the treat- 
ment of soft contributions and also higher-twist correc- 
tions appear to be significant. We estimate an irreducible 
theoretical uncertainty for this region (for the NLO cal- 
culation for a given DA) as ±5% but the accuracy of 
the duality approximation itself is difficult to quantify. 
It has to be tested. A sensitivity to the pion DA in this 
region is, on the other hand, mostly limited to the second- 
order coefficient ai in the Gegenbauer expansion which 
can be calculated rather precisely on the lattice. With 
this input, comparing the calculations of the 7r°7*7 tran- 
sition form factor with the data one will be able to get 
a much better understanding of the theoretical accuracy 
than presently available. This will be very important in 
a broader context, e.g. for QCD calculations of weak 
decays of B-mesons which involve similar relatively low 
scales. Especially with new data from the BES III col- 
laboration who plan a significant improvement over the 
Cleo, BaBar and Belle results in this region Q there can 
be quantitative progress. 



The existing three data points around Q 2 — 10 GeV 2 
have small error bars for both BaBar and BELLE, so that 
also the difference between the two experiments in this 
region is the most statistically significant. The large val- 
ues for g 2 F 7r o 7 . 7 (Q 2 ) ~ 0.019 in the Q 2 ~ 9 - 12 GeV 2 
region as compared to Q 2 F^o 1 * 1 {Q 2 ) ~ 0.015 at Q 2 ~ 
4—6 GeV 2 , reported by BaBar [3j , are difficult to de- 
scribe in the framework of QCD collincar factorization 
without invoking unconventional models of the pion DA 
with large end-point enhancements, see Ref. [l| for sev- 
eral examples and the discussion. The new data from 
BELLE Q 2 F 7T o r ~ f (Q 2 ) ~ 0.016 at Q 2 ~ 9 - 12 GeV 2 
are much easier to accommodate within a "standard" 
scenario: e.g the "holographic" model inspired by gauge- 
gravity correspondence provides a good fit, see above. 
This discrepancy must be clarified. Provided the Gegen- 
bauer coefficient ai is determined from lattice calcula- 
tions and/or fits to the data at lower Q 2 , experimental 
data at Q 2 ~ 10 GeV 2 will allow one to obtain quanti- 
tative information (constraints) on the shape of the pion 
DA beyond the value of its second moment. E.g. the 
coefficient 04 in the Gegenbauer expansion can be deter- 
mined in this way. 



3. large Q 2 > 20 GeV 2 



The transition form factor is expected to approach the 
scaling behavior in this region. To make this statement 
more quantitative, consider the following ratio: 



R 



(40/20) 



F^.^Q 2 = 40 GeV 2 ) 
'F n o rj (Q 2 = 20 GeV 2 ) 

We obtain (cf. Fig. QJ 

1.0K%o/20) <l-07, 



(2) 



(3) 



where the lower number corresponds to the BMS 
model [5| and the upper one to the "flat" DA. There 
is, on the other hand, no reason to expect that the form 
factor at Q 2 > 20 GeV 2 is close to the Brodsky-Lepage 
limit 



Q F„o~* 7 (Q ) 



2^Q 2 



= 0.185 



(4) 



because the asymptotic behavior is achieved in QCD very 
slowly. Thus, this value is much less constrained. For the 
range of the models shown in Fig. Q]we obtain e.g. 



o.i63 < [g 2 ^o 7 , 7 (g 2 )] 



5 2 =30 GeV 2 



< 0.245 (5) 



(all numbers in GeV) . Much of the excitement produced 
by the BaBar experimental data was due to an over- 
interpretation, from our point of view, of the power-law 
fit to the form factor Q 2 F n a , (Q 2 ) ~ 0.182(Q 2 /10) 025 
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provided in Ref . Q . This fit suggests a large scaling vio- 
lation i?f 4 a /2o) ~ 1-2 which is significantly outside of the 
range in Eq. ([3]) and difficult to reconcile with asymptotic 
freedom. The fit is, however, dominated by the data at 
lower momentum transfers so that such a rise at large 
Q 2 is not warranted. The same data can be described as 
(cf. [2]) 



Q'F^iQ 1 ) = AQ 2 /{B + Q Z 



(6) 



with the values of the parameters A = 0.230 ± 0.008 GeV 
and B = 2.57 ± 0.31 GeV 2 . Whereas the quality of this 
fit, x 2 — 1-73/d.o.f. is, admittedly, worse than of the 
BaBar power-law parametrization, \ 2 — 1.04/d.o.f., it 
shows that the evidence against the scaling behavior is 
at the level of less than two standard deviations. 

Since there are strong theoretical arguments that the 
scaling violation should not exceed a few percent, cf. 
Eq. ©, and since experimental errors are currently much 
larger, we prefer to average the existing data over the 
Q 2 > 20 GeV 2 region using kernel density estimation. 
Taking into account three data points with the largest 
values of Q 2 one obtains 



[Q 2 F,o 7 , 7 (g 2 )] 
[Q 2 F,o 7 , 7 (g 2 )] 



BaBar 

Q 2 >20 GoV 2 
BELLE 

Q 2 >20 GoV 2 



= 0.238 ± 0.029, 



0.204 ± 0.026 (7) 



for BaBar and BELLE experiments, respectively. We ob- 
serve that the two numbers are consistent to each other 
within 1.5a and are both somewhat larger than the QCD 
calculation with the asymptotic DA (the lower number in 
Eq. ([5])). This result supports the common wisdom that 
the pion DA at low scales is broader than the asymptotic 
distribution. A much higher accuracy is needed, how- 
ever, to discriminate between different models and pro- 
vide quantitative constraints. Future experimental data 
in the large-Q 2 region would also be very interesting for 
the comparison of time-like and space-like transition form 
factors which are predicted to be equal at asymptotically 
large momentum transfers. The corresponding discussion 
goes beyond the tasks of this short note. 

Finally, we want to make a few general remarks. First, 
one has to be very careful in comparing the statements 
on the shape of the pion DA that are obtained in dif- 
ferent theoretical approaches. For example, we find, in 
apparent contradiction to Ref. 0] , that a "flat" pion DA 
^flat _ ^ f a il s t describe the BaBar data, as it strongly 
overestimates the form factor at moderate Q 2 , see Fit* Q] 
The reason for this disagreement is that in Ref. [7[ a 
very large (nonstandard) nonperturbative correction is 
assumed that is beyond the operator product expansion. 
It is this assumption and not the shape of the pion DA 
that is crucial to obtain a strong scaling violation. 

Second, we shortly comment on the relation of our re- 
sults with those of the Bochum-Dubna group (for the 



latest update see which are obtained in the same ap- 
proach. One difference is that we do not include NNLO 
perturbative corrections induced by the running cou- 
pling and prefer to work consistently to the NLO ac- 
curacy. The reason is that terms corresponding the the 
running coupling (the so-called large-/?o approximation) 
usually strongly overestimate the full perturbative cor- 
rection. In addition there is danger of double count- 
ing with twist-four contributions because of QCD renor- 
malons, see [l(| for a review. Another difference is that 
we use somewhat larger values of the so-called Borel pa- 
rameter, M 2 ~ 1 — 2 GeV 2 vs. M 2 ~ 0.7 - 0.9 GeV 2 see 
Ref. [l[. Nevertheless, the numerical difference between 
our calculations (for a given DA) does not exceed 5-8% at 
Q 2 =4 — 8 GeV 2 , and appears to be negligible at larger 
momentum transfers. 

Note that in difference to Ref. [9| we do not have any 
preferred model of the pion DA. Our analysis is entirely 
devoted to the possibility of extraction of the parameters 
of the DA (Gegenbauer coefficients) from the experimen- 
tal data using collinear factorization, in full similarity 
with the determination of conventional parton distribu- 
tions from the data on inclusive reactions. 

We do not understand the claim made in Ref. [|| that 
parametrizations of the pion DA introduced in our work 
1] fail to describe the BaBar data: E.g. model II shown 
in Fig. [T] leads to \ 2 — 1-06 per data point, where the 
experimental errors and theoretical uncertainties of the 
calculation for a given DA are added in quadrature. For 
the central values of the parameters and not taking into 
account theoretical uncertainties we obtain for this model 
\ 2 = 1.34 which is also acceptable. 

To summarize, the measurements reported by the 
BELLE collaboration jH are by no means a game 
changer, but should somewhat take the heat off theo- 
rists struggling to invent a nonperturbative mechanism 
capable to postpone (or even invalidate, in the most dras- 
tic scenarios) the onset of QCD factorization in the pion 
transition form factor. The new generation of experimen- 
tal data that is expected to come from super-B factories 
will settle this question and allow the 7r°7*7 transition 
to serve its purpose as the gold-plated reaction in the 
theory of hard exclusive reactions. We have to warn, 
however, that this process alone is not sufficient for an 
unambiguous determination of the pion DA: one needs 
to do a global fit of all available hard exclusive processes 
including pions and eventually also include the lattice 
data. 
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